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ABSTRACT 

Nonrotating, zero metallicity stars with initial masses 140 260 Mq are expected to end their 

lives as pair-production supernovae (PPSNe), in which an electron-positron pair-production instability 
triggers explosive nuclear burning. Interest in such stars has been rekindled by recent theoretical studies 
that suggest primordial molecular clouds preferentially form stars with these masses. Since metal enrich- 
ment is a local process, the resulting PPSNe could occur over a broad range of redshifts, in pockets of 
metal-free gas. Using the implicit hydrodynamics code KEPLER, we have calculated a set of PPSN light 
curves that addresses the theoretical uncertainties and allows us to assess observational strategies for 
finding these objects at intermediate redshifts. The peak luminosities of typical PPSNe are only slightly 
greater than those of Type la, but they remain bright much longer (^^ 1 year) and have hydrogen lines. 
Ongoing supernova searches may soon be able to limit the contribution of these very massive stars to 
^ 1 % of the total star formation rate density out to z w 2, which already provides useful constraints for 
theoretical models. The planned Joint Dark Energy Mission satellite will be able to extend these limits 
out to z « 6. 

Subject headings: supernovae: general - galaxies: stellar content - stars:early-type 



1. INTRODUCTION 

Determining the cosmological evolution of the elements 
and elucidating the processes responsible for their creation 
are key aspects of understanding the history of the uni- 
verse. A wide spread in abundances is observed in stars, 
but the lowest metallicity galaxies studied so far are en- 
riched to substantial values of ~ 0.02.^© {e.g., Searle & 
Sargent 1972), and the diffuse medium in galaxy clusters is 
enriched to levels above 1/3 Zq {e.g., Renzini 1997). Simi- 
larly, quasar absorption line studies have uncovered heavy 
elements in even the lowest-density regions of the inter- 
galactic medium (IGM) {e.g., Schaye et al. 2003; Aracil 
et al. 2004) at levels that are roughly constant out to 
the highest redshifts probed (Songaila 2001; Pettini ct al. 
2003). And though direct observations of metal-poor ha- 
los stars have uncovered a number of unusual abundance 
patterns {e.g., Ryan, Norris, & Beers 1996; Christheb et 
al. 2002; Cayrel et al. 2004), no metal-free star has been 
detected to date. 

Thus a clear signature of the transition from a metal-free 
to an enriched universe has so far remained unobserved. 
Furthermore, there are good reasons to believe that this 
transition may have marked a drastic change in the charac- 
teristics of star formation. Recent theoretical studies have 
shown that collapsing primordial clouds fragmented into 
clumps with typical masses ~ 10^ M0 {e.g., Nakamura & 
Umemura 1998; Abel, Bryan, & Norman 2000; Bromm et 
al. 2001) and that accretion onto proto-stellar cores within 
these clouds was very efficient (Ripamonti et al. 2002; Tan 
& McKee 2004). The implication is that the initial mass 
function (IMF) of metal-free "Population III" stars may 



have been biased to higher masses than observed today. 

Since only a trace level of heavy elements will drasti- 
cally increase fragmentation, the formation of very massive 
stars (VMS) is unlikely beyond a "critical metallicity" of 
^ciit ^ 10~^ {e.g., Schneider et al. 2002; Bromm & Loeb 
2003). This means that the metals produced by the first 
few supernovae (SNe) were sufficient to halt very massive 
star formation, implying that this mode may have only 
operated at extremely high- redshifts z ^ 15 {e.g., Mackey, 
Bromm, & Hernquist 2003). Yet, such simple estimates ig- 
nore the substantial propagation times necessary for heavy 
elements to be mixed into pristine material. In fact, more 
detailed analyses show that it is almost impossible to 
raise the IGM metallicity above Zait in a homogeneous 
way (Scannapieco, Schneider, & Ferrara 2003, hereafter 
SSF03). Thus primordial metal enrichment was essentially 
a local process, which began in the densest "high-sigma" 
regions and gradually spread into to the more common 
lower-density regions as time progressed. This inhomo- 
geneous distribution parallels models of further IGM en- 
richment by second-generation objects {e.g., Scannapieco, 
Ferrara, & Madau 2002; Theuns et al. 2002; Fujita et al. 
2004; Cen, Nagamine & Ostriker 2004) and raises the pos- 
sibility that metal-free star formation may extend well into 
the observable range, in pockets of gas where metals have 
not yet propagated. 

Such primordial star-forming regions would naturally be 
isolated. Similarly, metal-free objects would tend to be 
just large enough for gravitational collapse to overcome 
the IGM thermal pressure, but small enough to have no 
direct progenitors above a cosmological Jeans mass. After 
reionization, such pristine star formation would be fun- 
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damcntally different tfian that talcing place in minihalos 
at very liigh redshifts, which depends sensitively on the 
presence of initial H2 {e.g., O'Shea et al. 2005). In small 
objects, molecular hydrogen is easily photodissociated by 
11.2-13.6 cV photons, to which the imi verse is otherwise 
transparent. Thus the emission from the first stars quickly 
destroyed all avenues for cooling by molecular line emis- 
sion (Dekel & Rees 1987; Haiman, Rees, & Loeb 1997; 
Ciardi, Ferrara, & Abel 2000). This quickly raised the 
minimum virial temperature necessary to cool effectively 
to approximately 10** K, although the precise value of this 
transition is the subject of debate {e.g., Glover & Brand 
2001 Yoshida et al. 2003) and is somewhat dependent on 
the level of the high-redshift X-ray background (Haiman, 
Rees, & Loeb 1996; Oh 2001; Machack, Bryan, & Abel 
2003). 

Later metal-free star formation is much more straight- 
forward, and is confined to halos with Tvir > 10^ K. In 
this case efficient atomic line cooling establishes a dense 
locally-stable disk, within which non-equilibrium free elec- 
trons catalyze H2. Unlike in less massive halos, H2 for- 
mation in these objects is largely impervious to feedback 
from external UV fields, due to the high densities achieved 
by atomic cooling (Oh & Haiman 2002). As discussed 
in SSF03, these lonely regions of late primordial star for- 
mation may have already been detected as a subgroup of 
high-redshift Lyman-alpha emitters, although distinguish- 
ing these objects from Population II stars is extremely 
difficult (Dawson et al. 2004). 

In this investigation we explore an alternative approach 
to search for such VMS: the direct detection of the re- 
sulting SNe. While many of the > 100 Mq stars 
would have died forming black holes, progenitors within 
the 140 — 260 Mq mass range end their lives in tremen- 
dously powerful pair- production supernovae (PPSNe). In 
these SNe, e+e" pair creation softens the equation state 
at the end of central carbon burning, leading to rapid col- 
lapse, followed by explosive burning up to the iron group 
{e.g., Ober, El Eid, & Fricke 1983; Bond, Arnett, Carr 
1984; Heger & Woosley 2002, hereafter HW02). The more 
massive the star, the higher the temperature at bounce 
and the heavier the elements that are produced by nuclear 
fusion. In all cases, the star is completely disrupted and 
both the ejection energies (~ 3 — 100 x 10^^ ergs) and 
masses (140 — 260 M©) are enormous, leading to observed 
properties qualitatively different from those of Type-la and 
typical core collapse SNe. In turn, these properties, such as 
bright UV colors, allow for efficient intermediate-redshift 
searches for very massive star formation. 

Recently, Wise & Abel (2005) studied the feasibility of 
detecting PPSNe at z ^ 20 with the James Webb Space 
Telescope, using extremely approximate lightcurves with 
luminosities that were taken to be proportional to the total 
^^Ni mass. A similar study was carried out by Weinmann 
& Lilly (2005), who also estimated number counts at z = 5 
using a single light curve that was taken from Heger ct 
al. (2002), which unfortunately contained a normalization 
error. ^ Limited studies of lightcurves from PPSNe have 
also been carried out by Woosley & Weaver (1982), who 
considered progenitors with only two masses and found 
features similar to those described below, and by Herzig 

®The curves in this conference proceeding are overluminous by a 



et al. (1990) who modeled the disruption of a 61 Mq Wolf- 
Rayet star, generating a relatively faint SN. Here we use 
the KEPLER code to directly construct a suite of PPSN 
models that addresses the range of theoretical possibilities 
and accounts for the detailed physical processes involved 
up until the point at which the SN becomes optically thin. 
These curves allow us to identify the properties that most 
distinguish PPSNe from Type la and core collapse SNe, 
and discuss their detectability at 2 ^ 6 in current and 
future surveys. 

The outline of this paper is as follows. In §2 we use 
the implicit hydrodynamical code KEPLER to construct 
a suite of detailed PPSN progenitor models. In §3 we gen- 
erate simulated spectra and lightcurves from these mod- 
els, highlighting the key features that distinguish PPSNe 
from other SNe. In §4 we apply these lightcurves to deter- 
mine the constraints on PPSNe obtainable by existing and 
planned SNe searches. In §5 we discuss the properties of 
the most likely environments of PPSNe, and we summarize 
our results in §6. 

2. MODELS OF PAIR-PRODUCTION SUPERNOVAE 

The progenitors of PPSNe are confined to a range of ini- 
tial masses that arc large enough to collapse via the pair- 
production instability, but small enough for this collapse 
to be reversed by nuclear burning. In the non-rotating 
case, these limits are clearly defined and well understood 
numerically (Ober, El Eid, & Fricke 1983; HW02; Heger 
et al. 2003). In stars above ^ 100 Mq, after central helium 
burning, the core reaches a sufficiently high temperature 
and entropy to produce e+e~ pairs, converting some of the 
internal energy of the gas into rest mass, which softens the 
equation of state to an adiabatic index below the critical 
value of 4/3 (Barkat, Rakavy, & Sack 1967; Bond, Arnett, 
& Carr 1984). The result is a runaway collapse that leads 
to explosive burning in the carbon-oxygen core. For stars 
with initial masses more than about 140 Mq, the energy 
released is sufficient to completely disrupt the star. 

A shock moves outward from the edge of the core, ini- 
tiating the supernova outburst when it reaches the stellar 
surface. Just above the 140 Mq limit, weak silicon burn- 
ing occurs and only trace amounts of radioactive '"'^Ni are 
produced. It is this ^^Ni that powers the late-time super- 
nova light curves. The amount of ^^Ni produced increases 
in larger progenitors, and in 260 M© progenitors up to 50 
Mq may be synthesized, ~ 100 times more than in a typ- 
ical Type la supernova. For stars with masses above 260 
Mq, however, the onset of photodisintegration in the cen- 
ter imposes an additional instability that collapses most 
of the star into a black hole (Bond, Arnett & Carr 1984; 
HW02; Fryer, Woosley, & Heger 2001). 

Pair-production supernovae are among most powerful 
thermonuclear explosions in the universe, with a total en- 
ergies ranging from 3 x 10^^ ergs for a 140 Mq star (64 
helium Mq core) to almost 100 x 10^^ ergs for a 260 Mq 
star (133 Mq helium core; HW02). We use the implicit 
hydrodynamical code KEPLER (Weaver, Zimmerman, & 
Woosley 1978) to model the entire evolution of the star 
and the resulting light curves. Since KEPLER only imple- 
ments gray diffusive radiation transport with approximate 
deposition of energy by gamma rays from radioactive de- 
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Table 1 

Properties of PPSN progenitor models 



Model 


Mho ( Mn) 


Mn ( Mn) 


Msewi ( Mn,) 


R (lO^'^cm) 


£^nn (lO'^^ergs) 


150-W 


70 


3.5(-4) 


4.2(-2) 


3.9 


6.9 


150-1 


46 


l.l(-4) 


6.3(-2) 


16 


9.2 


150-S 


49 


0.86 


8.6(-2) 


26 


8.5 


200-W 


97 


2.7(-6) 


3.3 


0.68 


29.5 


200-1 


58 


8.0(-6) 


5.1 


2.8 


36.5 


200-S 


89 


0.34 


2.2 


29 


29.1 


200-S2 


78 


4.75 


0.82 


20 


18.7 


250-W 


123 


3.1(-6) 


6.2 


0.58 


47.2 


250-1 


126 


9.1(-6) 


32 


4.0 


76.7 


250-S 


113 


1.34 


24.5 


26 


64.6 



cay of -"^^Ni and ^^Co (Eastman, Wooslcy, & Weaver 1993), 
the light curves obtained can only be followed as long as 
the SN is optically thick an there is a reasonably well- 
defined photosphere. Also during the very earliest stages 
of the supernova (shock break out), KEPLER underesti- 
mates the color temperature by about a factor 2 (Blinnikov 
et al. 2003). Finally, we employ a mixing procedure to sim- 
ulate, in 1-D, the effects of mixing by Rayleigh- Taylor in- 
stabilities after the supernova shock, which was calibrated 
using spectra of observed SN types. The main effect of 
this mixing for the light curves is the spreading of some 
^^Ni into the envelope. 

A supernova can be bright either because it makes a lot 
of radioactive "''''Ni (as in Type la supcrnovac) or because 
it has a large low density envelope and large radius (as in 
bright Type II supernovae). More radioactivity gives more 
energy at late times when the supernova has expanded 
enough to become transparent (R 10^^ cm). A larger 
initial radius means that the star cxpcric;nc;c;s less adiabatic; 
cooling as it expands to that radius, resulting in a higher 
luminosity at early times. Here, the most important fac- 
tors in determining the resulting light curves are the mass 
of the progenitor star and the efficiency of dredge-up of 
carbon from the core into the hydrogen envelope during or 
at the end of central helium burning. The specifics of the 
physical process encountered here are unique to primor- 
dial stars. Lacking initial metals, they have to produce 
the material for the CNO cycle themselves, through the 
synthesis of ^^C by the triple- alpha process. Just enough 
^^C is produced to initiate the CNO cycle and bring it 
into equilibrium: a mass fraction of 10~® when central hy- 
drogen burning starts, and a mass fraction ~ 10"'' during 
hydrogen-shell burning. 

At these low values the entropy in the hydrogen shell 
remains barely above that of the core, and the steep en- 
tropy gradient at the upper edge of the helium core that 
is typical for metal-enriched helium-burning stars is ab- 
sent. The high masses of the stars also imply that a large 
fraction of the pressure is due to radiation. This is well 
known to facilitate convection. For these reasons, the cen- 
tral convection zone during helium burning can get close, 
nip at, or even penetrate the hydrogen-rich layers. Once 
such mixing of high-temperature hydrogen and carbon oc- 
curs, the two components burn violently, and even without 



this rapid reaction, the hydrogen burning in the CNO cycle 
increases proportionately to the additional carbon. Thus 
mixing of material from the helium-burning core, which 
has a carbon abundance of order unity, is able to raise the 
energy generation rate in the hydrogen-burning shell by 
orders of magnitude over its intrinsic value. 

This mixing has two major effects on the PPSN progen- 
itor: first, it increases the opacity and energy generation 
in the envelope, leading to a red-giant structure for the 
presupernova star, in which the radius increases by over 
an order of magnitude. Second, it decreases the mass of 
the He core, consequently leading to a smaller mass of ^^Ni 
being synthesized and a smaller explosion energy. The for- 
mer effect increases the luminosity of the supernova, espe- 
cially at early times, while the latter effect can weaken 
it. Just how much mixing occurs is uncertain from stel- 
lar modeling at the present, though it has been studied 
by several authors {e.g., Heger, Woosley & Waters 2000; 
Marigo et al. 2001). In the present paper we account for 
these uncertainties by employing different values of con- 
vective overshooting - whose presence is clear, but whose 
quantitative interaction with the burning is unknown. 

A suite of representative models is chosen to address the 
expected range of presupernova models - from blue super- 
giant progenitors with little or no mixing to well-mixed 
red hypergiants, which can have pre-SN radii of 20 AU or 
more. These models arc summarized in Table 1, in which 
the names refer to the mass of the progenitor star (in units 
of M0) and the weak (W), intermediate (I), or strong (S) 
level of of convective overshoot. Note, that in all cases the 
progenitor star is nonrotating, and all models employ the 
same procedure to mix ^^Ni into the envelope after the 
supernova shock. 

The KEPLER code can be used to compute approx- 
imate light curves and has been validated both against 
much more complex and realistic codes such as EDDING- 
TON and observations of a prototypical Type II-P super- 
nova, SN 1969L (Weaver and Woosley 1980; Eastman et 
al 1994). Its deficiency is that it is a single temperature 
code using flux-limited radiative diffusion. It functions 
very well, however, in situations where the light curve is 
recombination dominated (as on the plateau of Type II- 
P supernovae), in the radioactive peaks of Type I super- 
nova, and in the tails of Type II-P supernovae. Gamma- 
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Fig. 1. — Luminosities (top row), effective temperatures (center row), and photospheric radii (bottom row) of PPSNe for ten different 
representative models. Each sets of panels are labeled by the total mass of the progenitor star, and models with weak (dashed), intermediate 
(solid) and strong (dot-dashed) convective overshoot are shown. Finally, the dotted lines in the central panels correspond to model 200-S2. 
See text and Table 1 for details. 



deposition has been calibrated against a Monte Carlo code 
(Pinto and Woosley 1988) and the code was successfully 
used to predict the behavior of SN 1987A (Woosley, Pinto, 
& Ensman 1988), though mixing was later added as an es- 
sential ingredient. Because thermal equilibrium and black 
body spectra are assumed, the bolometric luminosity is 
known much better than the color photometry. Thus B 
and V magnitudes should be approximately correct, but 
U magnitudes are much more uncertain. 

The resulting luminosities, effective temperatures, pho- 
tospheric radii are shown in Figure 1. As the shock moves 
toward the low-density stellar surface, its energy is de- 
posited into progressively smaller amounts of matter. This 
results in high velocities and temperatures when the shock 
reaches the stellar surface, causing a pulse of ultraviolet 
radiation with a characteristic timescale of a few minutes. 
This "breakout" phase is by far the most luminous and 
bluest phase of the PISN burst, but its very short dura- 
tion makes it difficult to use in observational searches. In 
fact, the analog of this phase in conventional SN has so 
far only been indirectly detected in SN 1987A (Nayozhin 
1994; Hamuy et al. 1988; Catchpole et al. 1988). 

Following breakout, the star expands with i?phot ini- 
tially proportional to time. Though a small fraction of 



the outer mass may move much faster, the characteristic 
velocity of the photosphere during this phase is a modest 
V = {2KE/Mf/'^ ~ (lO53gj.gs/2OOM0)i/2 _ 500O km/s, 
because of the very large mass participating in the ex- 
plosion. During the expansion, the radiation-dominated 
ejecta cool adiabatically, with T approximately propor- 
tional to , with an additional energy input from the 
decay of '^^Ni (if a significant mass was synthesized dur- 
ing the explosion) and hydrogen recombinations (when 
T ~ 10'*K). As the scale radius for this cooling is the ra- 
dius of the progenitor, the temperatures and luminosities 
are substantially larger throughout this phase in the cases 
with the strongest mixing. 

After ~ 50 days, the energy input from ^^Co decay be- 
comes larger than the remaining thermal energy (the ini- 
tial thermal energy deposited by the shock is mostly eaten 
away by the adiabatic expansion) . At this time the energy 
deposited by ^^Co in deeper layers that were enriched in 
^^Ni can diffuse out and these layers also become gradually 
more exposed as the outer parts of the supernova ejecta 
reconibine and become optically thin. For stars that were 
compact to begin with, this can cause a delayed rise to 
the peak of the light curve. For stars with larger radii, the 
radioactivity just makes a bright tail following the long 



^ 26 
25 



^ 29 



L PPSN(200-W) JL PPSN(200-I) JL PPSN(200-S) J 




28 
27 
26 

25 L PPSN(250-W) _|_ PPSN(250-I) _|_ PPSN(250-S) J 



2000 5000 



8000 5000 8000 

wave length [A] 



5000 8000 



Fig. 2. — Comparison of PPSNe spectra with those of Type la and Type II-P supernovae at 1 day (dashed), 10 day (dot-dashed), 100 days 
(sohd), and 200 days (dotted) since explosion. The Type la spectra were taken from observations of SN 1994D, the Type II-P spectra were 
taken from observations of SN 1999em, and the PPSNe spectra were computed from black body models with Tj-oior = T^fi- While the overall 
shape of our simplified PPSNe spectra are similar to those of SNe Type la and Type II, their late-time evolution is markedly different, as 
seem most clearly in the 200 day curves. 



plateau in emission from the expanding envelope. Even- 
tually, even the slow-moving inner layers recombine and 
there is no longer a well-defined photosphere. At this time 
the assumption of local thermodynamic equilibrium (LTE) 
breaks down, and more detailed radiative transfer calcu- 
lations are required, which are beyond the scope of our 
modeling here. The SN is fainter and redder during this 
phase, however, and thus difficult to detect at cosmological 
distances in optical and near infrared (NIR) surveys. 

3. SPECTRAL PROPERTIES AND LIGHTCURVES 

From the luminosities and effective temperatures in Fig- 
ure 1, we calculated approximate PPSNe spectra, assum- 
ing a black body distribution with the color temperature 
equal to the effective temperature, that is the specific lu- 
minosity 

L.it)=L{t) (1) 

where x = v/vq and vq = 2.1 x fO^^T4.cff Hz, with 
r4_off = Toff/10'* K. Recall that the peak frequency in this 
case occurs at 5.9 x 10*^r4_ofT Hz, which corresponds to 



a wavelength of 5100TjJjj A. The resulting spectra are 
plotted at representative times of 1, 10, 100, and 200 days 
in Figure 2, in which we also include comparison spectra 
from observations of Type la and Type II (core collapse) 
SNe. In particular, the Type la curves were taken from 
observations of SN 1994D by Patat et al. (1996) and Fil- 
ippenko (1997) as compiled by Filippenko (1999), at 1, 10 
and 128 days, and 233 days since explosion (assumed to be 
13 days before B-max) normalized to Patat et al. (1996) 
and Cappelaro et al. (1997) V-band light curves at 1, 10, 
100, and 200 days, with an assumed distance of 13.7 Mpc 
(see also Richmond et al. 1995). The Type II curves, on 
the other hand, were taken from spectra of the Type II-P 
SN 1999em by Elmhamdi et al. (2003) at 9, 10, 113, and 
168 days since explosion, again normalized to their V-band 
light curves at 1, 10, 100, and 200 days, with an assumed 
distance of 7.8 Mpc. 

The most striking feature from this comparison is that 
despite enormous kinetic energies of ~ 50 x 10^^ ergs, the 
peak optical luminosities of PPSNe are similar to those of 
other SNe, even falling below the la and II curves in many 
cases. This is because the higher ejecta mass produces 
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Fig. 3. — Comparison of light curves of a SN Type la, a SN Type II-P, a bright SN Type II-L, and PPSNe models with varying progenitor 
masses. In all cases the solid lines are absolute V-band AB magnitudes, the dot-dashed lines are the absolute B-band AB magnitudes, and 
the dashed lines are the absolute U-band AB magnitudes. 




a large optical depth and most of the internal energy of 
the gas is converted into kinetic energy by adiabatic ex- 
pansion. Furthermore, the overall spectral shapes of the 
PPSN curves are not unlike those of more usual cases. In 
fact, pair-production supernovae spend most of their lives 
in the same temperature range as other SNe and therefore 
exhibit a similar range in colors throughout their evolu- 
tion. Lastly, many of line features in the observed spectra 
would also be present in more detailed models of PPSNe. 
In particular, because of their hydrogen envelopes, PPSN 
spectra should contain hydrogen lines similar to those in 
Type-II SNe. Clearly, then, PPSNe will not be obviously 
distinguishable from their more usual counterparts "at 
first glance." 

A closer comparison between spectra, however, uncov- 
ers two key features that are uniquely characteristic to 
PPSNe. The first of these is a dramatically extended in- 
trinsic decay time, which is especially noticeable in the 



models with the strongest enrichment of CNO in the en- 
velope. This is due to the long adiabatic cooling times 
of supergiant progenitors, whose radii are ~ 20 AU, but 
whose expansion velocities are similar or even less than 
those of other SNe. Second, PPSNe are the only objects 
that show an extremely late rise at times > 100 days. This 
is due to energy released by the decay of ^^Co, which un- 
like in the Type la case, takes months to dominate over the 
internal energy imparted by the initial shock. In this case 
the feature is strongest in models with the least mixing 
and envelope enrichment during helium burning, as these 
have the largest helium cores and consequently the largest 
^^Ni masses. Note, however, that neither of these features 
is generically present in all PPSNe, and both can be ab- 
sent in smaller VMS that fail expand to large sizes through 
dredge-up and do not synthesize appreciable amounts of 
^^Ni. In the 150-W case, for example, the luminosity de- 
cays monotonically on a relatively short time scale, pro- 
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ducing spectra not dissimilar to the comparison Type II 
curves from SN 1999em. In fact this 150 Mq SN shares 
many similarities with its smaller-mass cousin: both are 
SNe from progenitors with radii ~ 10^^ cm and in both 
^^Ni plays a negligible role. 

In Figure 3 we examine the temporal evolution of PP- 
SNe in more detail by plotting absolute AB light curves 
of our models at three representative wavelengths: 5500 
A, corresponding to the central wavelength of the V-band; 
4400 A, corresponding to the B-band; and 3650 A, corre- 
sponding to the U-band. We focus on blue wavelengths as 
it is features in these bands that will be redshifted into the 
optical and NIR at cosmological distances. Again, for com- 
parison, we also include in this figure observed light curves 
for SN Type la and Type II. In the Type la case the curves 
are again taken from observations of 1994D by Patat et al. 
(1996) and are supplemented at late times by data from 
Cappelaro et al. (1997). In the Type II case, we consider 
both Type II-P and Type II-L SNe. Our Type II-P curves 
are taken from observations of 1999em by Elmhamdi et 
al. (2003), and our Type II-L curves are taken from obser- 
vations of the very bright supernova 1979C, as compiled 
in de Vaucouleurs et al. (1981) and Barbon et al. (1982), 
with an assmned distance of 17.2 Mpc (Freedman et al. 
1994). 

From this point of view, the characteristic features of 
PPSNe are even more striking, and the relation with pro- 
genitor structure is clear. Compact models with weak 
dredge-up are dominated by a late-time ^^Ni bump; ex- 
tended hypergiant progenitor models with strong mixing 
show a significant blue early phase during the expansion 
followed by a longer, redder phase in which the luminosity 
is roughly constant for ^ 200 days; and intermediate mod- 
els are a hybrid of the two. Finally, in cases such as the 
compact 150-W blue supergiant model, in which dredge- up 
does not occur and little ®^Ni is synthesized, both of these 
phases are absent, leaving an object remarkably reminis- 
cent of more typical core collapse SNe, as exemplified by 
the Type II-P lightcurves. Similarly, models such as 150-1 
and 150-S which also synthesize very little ^^Ni, but are 
preceded by larger progenitors, are easily confused with 
bright Type II-L SNe, particularly at early times. In no 
case, however, do PPSNe look anything like SNe Type la. 
In particular none of the pair-production models display 
the long exponential decay seen in the Type la curves, 
and all PPSNe contain hydrogen lines, arising from their 
substantial envelopes. 

4. PAIR-PRODUCTION SUPERNOVAE IN COSMOLOGICAL 

SURVEYS 

From the models developed in §3, it is relatively straight- 
forward to relate the star formation history of VMS to the 
resulting number of observable pair-production supernova. 
In this section and below we adopt cosmological param- 
eters of /i = 0.7, = 0.3, Vtf, = 0.7, and = 0.045, 
where h is the Hubble constant in units of 100 km s"-"^ and 
^rnj ^A, and fib are the total matter, vacuum, and bay- 
onic densities in units of the critical density {e.g., Spergel 
et al. 2003). 

Here we focus on three PPSN light curves, which bracket 
the range of possibilities: the faintest of all our models, 
150-W, in which both significant dredge-up and ^^Ni pro- 
duction are absent; an intermediate model, 200-1, in which 



little dredge-up occurred, but 5.1 Mq of ^^Ni were formed; 
and the model with the brightest lightcurves, 250-S, in 
which substantial dredge-up leads to an enormous initial 
radius of over 20 A.U., and the production of 24.5 Mq 
of ^^Ni causes an extended late-time period of high lumi- 
nosity. Accounting for redshifting, time dilation, and the 
appropriate luminosity and bandwidth factors, the specific 
flux for each of these models at a redshift z observed at a 
wavelength A and a time t after breakout (in the frame of 
the observer) is 
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where r{z) is the comoving distance associated with the 
supernova redshift. Note that this expression does not ad- 
dress the possibility of extinction by dust, which amounts 
to assuming that pristine regions remain dust-free thought 
the lifetime of the very massive PPSN progenitors stars. 
For any given PPSN model, we can then use eq. (2) to cal- 
culate t{\, -F™™, z) the total time the observed flux at the 
wavelength A from a SN at the redshift z is greater than the 
magnitude limit associated with the specific flux i^™'". Fi- 
nally, the total number of pair-production SNe shining at 
any given time with fluxes above F™'", per square degree 
per unit redshift is given by the product of the volume ele- 
ment, the (time-dilated) PPSN rate density, and the time 
a given PPSN is visible, that is 
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where the rate density p^{z) is the number of PPSNe per 
unit time per comoving volume as a function of redshift. 

As p^{z) is completely unknown, we adopt here two 
simple models. In the first model, we assume that 
metal-free star formation occurs at a constant rate den- 
sity, which we take to be 0.001 MQyr"-"^ Mpc~^, indepen- 
dent of redshift. In the second case, we assume that 
at all redshifts metal- free stars form at 1% of the ob- 
served total star formation rate density, which we model 
as logiolpf ^(2)/M0yr-iMpc-3] = _2.1 + 3.9 logio(l + 
2;) — 3.0 [logjQ(l + z)]^, a simple fit to the most recent mea- 
surements (Giavalisco et al. 2004; Bouwens et al. 2004). 
Finally, for both star formation models we assume that 1 
pair-production SN occurs per 1000 solar masses of metal 
free stars. This is consistent with typical estimates given 
in SSF03, in which N^, defined as the number of PPSNe 
per solar mass of stars formed, was found to vary between 
~ 10"'* and ~ 5 X 10~^ for a wide range of possible metal- 
free IMFs. 

Note that while our two simple PPSN rate densities 
are chosen such that they can be easily rescaled by the 
reader, they are nevertheless consistent with the range of 
z ^ 1 values predicted in more sophisticated models, as 
discussed in §5. Furthermore, the total amount of metals 
produced in our simple models is consistent with the ele- 
ment abundances observed in extremely metal-poor Galac- 
tic halo stars. Assuming a typical value of 200 Mq of met- 
als ejected per PPSN and integrating down to z = 2.5, one 
obtains values ~ 6 x 10^ Mq Mpc~^ for both of our PPSN 
rate density models. This corresponds to a mass fraction 
of baryons that have been processed by VMS of ^ 6 x 10~^, 
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Fig. 4. — Number of PPSNe per square degree per unit redshift above a given I-band magnitude, assuming 0.001 M0 yr-i Mpc-^ (solid 
lines) or 1% of the observed star formation rate density (dashed lines). The Iab = 26 cut taken in the upper rows approximately corresponds 
to the magnitude limit of the Institute for Astronomy Deep Survey which covered 2.5 deg^ (as shown by the dotted lines) from September 
2001 to April 2002. The Ij^g = 27 limit in the bottom panel corresponds to that of the ongoing COSMOS survey, which will survey an area 
of 2 deg^ (again indicated by the dotted lines). Note however, that the COSMOS survey itself is primarily focused on large-scale structure 
issues and will not be able to find PPSNe, as each pointing is visited only once. 



which is consistent with the 3 — 7 x lO^'^ hmits inferred 
by Oh et al. (2001) to explain the relative abundances in 
extremely metal-poor Galactic stars. 

The resulting observed PPSNe counts for these models 
are given in Figure 4 for two limiting magnitudes. In the 
upper panels, we take a Iab = 26 magnitude limit, appro- 
priate for the Institute for Astronomy (IfA) Deep Survey 
(Barris et al. 2004) , a ground-based survey that covered a 
total of 2.5 deg2 from September 2001 to April 2002. As we 
are interested in rare objects, this type of survey is more 
constraining that a more detailed, smaller-area searches 
such as the Hubble Higher z Supernova Search (Riess et 
al. 2004; Strolger et al. 2004). 

From this figure we see that existing data sets, if prop- 
erly analyzed, are easily able to place useful constraints on 
VMS formation at low redshifts. Given a typical PPSN 
model like 200-1 for example, the already realized IfA sur- 
vey can be used to place a constraint of < 1% of the total 
star formation rate density out to a redshift ~ 1. Simi- 
larly, extreme models such as 250-S can be probed out to 
redshifts ~ 2, all within the context of a recent SN search 
driven by completely different science goals. Note how- 
ever that these limits are strongly dependent on significant 
mixing in the SN progenitor or the production of ^^Ni, and 
thus models such as 150-W remain largely unconstrained 
by the IfA survey. 

In the bottom panels of Figure 4 we consider a limit- 
ing magnitude of Iab — 27, appropriate for the COSMOS 

'^see http:/ /www. astro. caltech.cdu/'^cosmos/ 
*see http://snap.lbl.gov/ 
®see http://destiny.asu.edu/ 



survey^, an ongoing project that will cover 2 deg^ using 
the Advanced Camera for Surveys on HST. Raising the 
limiting magnitude from Iab — 26 to Iab ~ 27 has the 
primary effect of extending the sensitivity out to slightly 
higher redshifts. This pushes the probed range from z < I 
to z < 1.5 in the 200-1 case and from z < 2 to z 3 in 
the 250-S case. Again this is all in the context of an on- 
going survey. Even with this fainter limiting magnitude, 
however, low-luminosity PPSNe like 150-W are extremely 
difficult to find, and remain largely unconstrained. 

This shortcoming is easily overcome by moving to NIR 
wavelengths. In Figure 5, we calculate the PPSN con- 
straints that would be obtained from three possible realiza- 
tions of the planned space-based Joint Dark Energy Mis- 
sion (JDEM). In the upper and central panels of this plot 
we adopt two limiting J-band magnitudes, appropriate for 
the two surveys that would be carried out with the Super- 
nova Acceleration Cosmology Probe (SNAP)^ realization of 
JDEM. In the upper panel we adopt a limit of Jab = 27.7, 
corresponding to the depth of the planned (700 deg^) wide- 
field SNAP survey. With 2 months of data, the alterna- 
tive Dark Energy Space Telescope (Destiny)^ realization 
of JDEM would cover 7.5 deg^ of sky with spectroscopic 
observations down to a similar magnitude limit. 

In these NIR surveys, dN^cg^ /dz is dramatically in- 
creased with respect to ground-based searches. This is 
due to the fact that for the majority of their lifetimes, 
the effective temperatures of PPSNe are just above the 
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Fig. 5. — Number of PPSNe per square degree per unit redshift above a fixed J or K-band magnitude, with limits appropriate for three 
possible realizations of the Joint Dark Energy Mission. Curves are as in Figure 4. Near infrared coverage allows for detection of PPSNe even 
after the lightcurves have been significantly redshifted. Top: The Jab = 27.7 limit taken in these panels corresponds to that of a single scan 
in the planned deep-field SNAP survey, which would cover 7.5 deg-^ (indicated by the upper dotted lines). A similar magnitude limit with the 
same coverage would be obtained with 2 months of grism data from the Destiny mission. The full planned 700 deg^ wide-field SNAP survey 
(lower dotted lines) also has the same limiting magnitude. Center: The Jab = 30.3 limit taken in these panels is that of the deep-field SNAP 
survey, which would be able to place constraints on PPSNe out to 2 > 5. Bottom: Curves corresponding to the JEDI mission, whose redder 

photometry coverage allows for detections down to = 28.7 over 24 deg^ with 2 months of data. This would place limits on PPSNe out 

to z ~ 6 and beyond. 



^ IQ^'^K recombination temperature of hydrogen, which 
corresponds to a peak black-body wavelength ^ 8000 A. 
This means that for all but the lowest redshifts, the ma- 
jority of the emitted light is shifted substantially redward 
of the I-band, which is centered at 9000 A. Thus moving 
to the J-band, which is centered at 12500 A, represents 
an exponential increase in the observed flux and allows for 
detections at significantly higher redshifts. 

In the 200-1 case, this results in number densities that 
are almost an order of magnitude higher than in the I- 
band. Combined with the larger area surveyed, this means 
that the wide-area SNAP survey could place constraints on 
this model that are roughly 300 times more stringent than 
those from the IfA survey, and these constraints would ex- 

^'^see http://jedi.nhn.ou.edu/ 



tend to z ^ 4. Similarly tight limits can be placed on the 
250-S model, but this time extending out to a redshift of 
6. In fact, even very faint PPSNe like 150-W would be 
sensitively probed out to z ^ 2.5. In the central panels 
we adopt a limit of Jab = 30, corresponding to the depth 
of the planned (7.5 deg^) SNAP deep-field survey. In this 
case, both the 200-1 and 250-S models would be well stud- 
ied out to z ^ 5, and the 150-W model would be easily 
detectable out to z ^ 4 despite its overall low luminosity 
and relatively short lifetime. 

Finally, in the lower panel, we consider an even redder 
24 deg^ survey with a limiting magnitude of Kab — 28.7, 
as appropriate for two months of observations from the 
Joint Efficient Dark-energy Investigation (JEDI)^'^ real- 
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ization of JDEM. As the K band is centered at 22000 A, 
such observations naturally push to even higher redshifts. 
Thus while reaching a limiting AB magnitude only slightly 
higher than the 700 deg^ SNAP survey, such a 24 deg^ is 
able to place the most stringent of any of the surveys con- 
sidered: constraining the 150-W model out to z ~ 4 and 
the more luminous models to z ^ 6 and beyond. 



5. THE ENVIRONMENTS OF PAIR-PRODUCTION 
SUPERNOVAE 

As cosmological enrichment is essentially a local pro- 
cess, certain environments are naturally more favorable for 
metal- free star formation. In SSF03, we showed that the 
transition from mctal-frcc to Population II stars was heav- 
ily dependent on the efficiency with which metals where 
mixed into the intergalactic medium. This efficiency de- 
pended in turn on the energy input into galactic outflows 
powered by PPSNe, which was parameterized by the "en- 
ergy input per unit primordial gas mass" ^^g^^, defined as 
the product of the fraction of gas in each primordial object 
that is converted into stars (//^^), the number of PPSNe 
per unit mass of metal- free stars formed (7V^ ) , the average 
kinetic energy per pair- production supernova (fkin), and 
the fraction of the total kinetic energy channeled into the 
resulting galaxy outflow (/wind)- 

Incorporating such outflows into a detailed analytical 
model of structure formation leads to the approximate re- 
lation that, by mass, the fraction of the total star forma- 
tion in metal-free stars at z = 4 is 
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Fig. 6. — Number of PPSNe per square degree per unit redshift per 
year for a wide range of models. As in Figs. 4 and 5, the solid 
and dashed curves assume Pop III star formation rate densities of 
0.001 Mq yr~^ Mpc"'' and 1 % of the observed star formation 
rate density, respectively. The shaded region covers the range of 
metal-free star formation rate density models considered in SSF03, 



with the weakest feedback model (£ 



III 



10 ) defining the upper 



end, and the strongest feedback model = 10 



2.5', 



.,g — , defining the 

lower end. An extrapolation of these star formation rate densities 
to 2 = leads to the crosshatched region. In all SSF03 models the 
highest rates occur at redshifts < 10. Finally, the starred point is 
the 2 = 15 estimate by Weinmann & Lilly (2005). 



where, as above, E is in units of 10^"'^ ergs per Mq of 
gas (see Figure 3 of SSF03 for details). Extrapolating the 
results in SSF03 to z = gives 

Fl"{z = 0) - 10-^-5(£f ^)-\ (5) 

These fractions can be related to the underlying popu- 
lation of stars by adopting fiducial values of fl^^ — 0.1 
for the star formation efficiency, which is consistent with 
the observed star formation rate density at intermediate 
and high redshifts (Scannapieco, Ferrara, & Madau 2002); 
fw = 0.3 for the wind efficiency, which is consistent with 
the dwarf galaxy outflow simulations of Mori, Ferrara, & 
Madau (2002); and iV± = 0.001, for the number of PPSNe 
per unit solar mass of primordial stars formed, which is the 
value assumed in §4. This gives Fl^^ values of 0.3(i?kin)^^ 
at z = 4 and 0.1(i?kin)^^ at z = 0, respectively. Or, in 
other words, for typical energies of 30 x 10^^ ergs per PP- 
SNe, ~ 1% of the star formation aX z — A and ~ 0.3% of 
the star formation at z = by mass could be in metal-free. 

In Figure 6 we show estimates of the number of SNe 
per deg^ per dz per year over the wide range of models 
considered in SSF03, extrapolating to z = 0. In all cases 
we assume that 1 PPSN forms per 1000 solar masses of 
metal-free stars, and for comparison we show the simple 
low-redshift estimates taken in the previous section. While 
the values in these models are uncertain, they nevertheless 



serve to illustrate two important points. First, even at the 
low redshifts probed by current surveys, the simple range 
of VMS formation rates considered in §4 lie well within 
the theoretically interesting range. Indeed for the weakest 
feedback cases several z ~ 1 PPSNe should have already 
been seen within the IfA Deep survey. Secondly, due to 
the decrease in dr/dz as a function of redshift, as well as 
the time dilation effect, the peak of dN^cgl dzdt is not in 
the z ^ 10 range targeted by the James Webb Space Tele- 
scope^ but rather at more moderate redshifts. This means 
that the constraints obtained from JDEM-type surveys are 
likely to reach comparable limits to higher-rcdshift efforts. 
In fact, it is extremely unlikely that PPSNe will be found 
at higher redshifts without similar detections at z < 6. 

Note that for the full range of models in Figure 6, metal- 
free star formation naturally occurs in the smallest galax- 
ies, just large enough to overcome the thermal pressure of 
the ionized IGM, but small enough not to be clustered 
near areas of previous star formation (SSF03). In our 
adopted cosmology, for an IGM temperature of 2 x lO'', 
the minimum virial mass is 8.4 x 10^(1 + z)~^^^ M© with 
a corresponding gas mass of 1.4 X 10^(1 + z)-3/2 ]VJg_ Tj^is 
means the total stellar mass of primordial objects is likely 
to be around ~ 10*Mq(1 + z)"^/^, many orders of 
magnitude below L^, galaxies. Thus in general blank-field 
surveys should be the best method for searching for PP- 
SNe, as catalogs of likely host galaxies would be extremely 
difficult to construct. 

Nevertheless, as VMS shine so brightly, a direct search 
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for primordial host galaxies is not a hopeless endeavor. 
In particularly the lack of dust in these objects and the 
large number of ionizing photons from massive metal-free 
stars leads naturally to a greatly enhanced Lyman alpha 
luminosity. Following SSF03 this can be estimated as 

L^=CL{l-fesc)QiH)M^, (6) 

where cl = 1.04 x 10~^^ ergs, /esc is the escape fraction 

of ionizing photons from the galaxy, which is likely to be 
^ 0.2 (see Ciardi, Bianchi, & Ferrara 2002 and references 
therein), and the ionizing photon rate Q{H) can be esti- 
mated as « 10''^ s^i Mq~^ (Schaerer 2002). This gives a 
value of La ^ 10^^(1 + z)^'^/^ergs s~^ which, if observed 
in a typical 1000 A wide broad band corresponds to an 
absolute AB mag ~ —23 + 3.81og(l + z), much brighter 
than the PPSNe themselves. However this flux would be 
spread out over many pixels and be more difficult to ob- 
serve against the sky than the point-like PPSNe emission. 
For further details on the detectability of metal- free stars 
though Lyman-alpha observations, the reader is referred 
to SSF03. 

Finally we note that the presence of pair-production su- 
pernovae is generic to stars that die within a given range 
of helium core masses (HW02), and does not directly de- 
pend on metallicity. However there are three important 
theoretical reasons that make metal-free stars the favored 
progenitors of PPSNe. First, as H2 becomes an ineffi- 
cient coolant below a typical density and temperature of 
~ 10^ cm~^ and ~ 100 K, fragmentation to masses smaller 
than ^ 10^ Mq is highly suppressed in gas of primordial 
composition {e.g., Abel, Bryan, & Norman 2002). Sec- 
ondly, metal-rich VMS are prone to opacity-driven radial 
pulsations, which are suppressed in stars with metallici- 
ties below ~ IQ-^Zq (Baraffe, Heger, & Woosley 2001). 
Lastly, the lack of metals greatly reduces the line-driven 
wind mass loss, which is predicted to decline with metal- 
licity as or faster (Kudritzki 2000; Vink et al. 2001; 
Kudritzki 2002). 

From an observational point of view, Figer (2005) car- 
ried out a detailed study of the IMF in the Arches cluster, 
which is large {M^ > 10'^ Mq), young (r = 2.0 - 2.5 Myrs), 
and at a well-determined distance, making it ideal for such 
studies. No stars more massive than 130 M0 were found 
in this Z ~ Zq cluster, although more than 18 were ex- 
pected. A similar ~ 150 M0 limit was found in the lower 
metallicity cluster R136 in the Large Magellanic Cloud 
(Weidener & Kroupa 2003), and while mass estimates of 
a few stars exceed this limit, these values are quite uncer- 
tain. The Pistol star, for example, may in fact be a tight 
binary or have recently experienced a merger with another 
star (Figer & Kim 2002). 

Intriguingly, the Arches cluster not only exhibits a cut- 
off at the highest masses but a flattened IMF above 50 M0 
(r = —0.9 rather than the —1.35 Salpeter slope). Could 
some of these stars have been quickly whittled down from 
larger-mass progenitors? The answer is unclear. What is 
clear, however, is that there is little observational or the- 
oretical evidence that the progenitors of PPSNe might be 
found in enriched environments. 

6. SUMMARY 



Astronomers naturally associate metal-free star forma- 
tion with extremely high redshifts. While the c^arly uni- 
verse contained no elements heavier than lithium, today 
stellar nucleosynthetic products are found in all measured 
Galactic halo stars, all nearby galaxies, and even in the 
low-density IGM. Yet it would be a mistake to conclude 
that such observations exclude metal-free star formation 
at moderate redshifts. Metal-enrichment is an intrinsi- 
cally local process that proceeds over an extended redshift 
range, and at each redshift, the pockets of metal-frcx- star 
formation are naturally confined to the lowest-mass galax- 
ies, which are small enough not to be clustered near ar- 
eas of previous star formation. As such faint galaxies are 
difficult to detect and even more difficult to confirm as 
metal-free, the hosts of PPSNe could easily be lurking at 
the limits of present-day galaxy surveys. 

Using the implicit hydrodynamical code KEPLER we 
have constructed a suite of lightcurves that address the 
theoretical uncertainties involved in modeling PPSNe. 
Here the most important factors are the mass of the pro- 
genitor star and the efficiency of dredge-up of carbon from 
the core into the envelope. In general, increasing the mass 
leads to greater ^^Ni production, which boots the late time 
SN luminosity. Mixing, on the other hand, has two major 
effects: it increases the opacity in the envelope, leading to 
a red giant phase that increases the early-time SN luminos- 
ity; and it decreases the mass of the He core, consequently 
leading to a somewhat smaller mass of ^^Ni being synthe- 
sized. Despite these uncertainties, PPSNe in general can 
be characterized by three key features: (1) peak magni- 
tudes that are brighter than Type II SNe and comparable 
or slightly brighter than typical SNe Type la; (2) very long 
decay times ^ 1 year, which result from the large initial 
radii and large masses of material involved in the explo- 
sion; and (3) the presence of hydrogen lines, which are 
caused by the outer envelope. Note, however, only this 
last feature is present in all cases, and in fact, the lowest 
mass PPSN models we constructed have lightcurves that 
are remarkably similar to those of SN Type II. 

Accounting for redshifting, time dilation, and appropri- 
ate luminosity factors we used these lightcurves to relate 
the overall very massive star formation rate density to the 
number of PPSNe detectable in current and planned su- 
pernova searches. Here the long lifetimes help to keep a 
substantial number of PPSNe visible at any given time, 
meaning that ongoing SN searches should be able to limit 
the contribution of VMS to 1 % of the total star forma- 
tion rate density out to a redshift of 2, unless both mixing 
and ^^Ni production are absent for all PPSNe. Such con- 
straints already place meaningful limits on the cosmologi- 
cal propagation of metals. 

The impact of future NIR searches is even more promis- 
ing, as the majority of the PPSN light is emitted at rest- 
frame wavelengths longward of ~ 8000 A. Thus planned 
NIR satellite missions such as JD EM would be over two or- 
ders of magnitudes more sensitive to PPSNe than present 
optical surveys, and able to probe redshifts beyond z ^ 6. 
In this case, even the dimmest PPSNe would be detectable 
out to z « 4. 

Although the peak of the metal-free star formation den- 
sity almost certainly occurred at extremely early times, 
there is much to be learned from PPSN searches at more 
moderate redshifts. In fact, due to volume and time- 
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dilation effects, the peak in the number of PPSNe per deg^ 
per dz per year is hkely to he well below z = 10. Further- 
more, the data sets necessary for such analyses are already 
being planned for and collected. While the properties PP- 
SNe are diverse, a singular conclusion can be drawn from 
our modeling. Searches for pair-production SNe at 2; ^ 6 
will dramatically increase our understanding of the history 
of cosmic enrichment, the nature of metal-free stars, and 
the evolution of gaseous matter in the universe. 
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